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Abstract 

We present a systematic theoretical study of the five smallest oligoacenes (naph- 
thalene, anthracene, tetracene, pentacene, and hexacene) in their anionic, neutral, 
cationic, and dicationic charge states. We used density functional theory (DFT) to 
obtain the ground-state optimised geometries, and time-dependent DFT (TD-DFT) 
to evaluate the electronic absorption spectra. Total-energy differences enabled us to 
evaluate the electron affinities and first and second ionisation energies, the quasipar- 
ticle correction to the HOMO-LUMO energy gap and an estimate of the excitonic ef- 
fects in the neutral molecules. Electronic absorption spectra have been computed by 
combining two different implementations of TD-DFT: the frequency-space method 
to study general trends as a function of charge-state and molecular size for the 
lowest-lying in-plane long-polarised and short-polarised vr — > vr* electronic transi- 
tions, and the real-time propagation scheme to obtain the whole photo-absorption 
cross-section up to the far-UV. Doubly-ionised PAHs are found to display strong 
electronic transitions of vr — > tt* character in the near-IR, visible, and near-UV spec- 
tral ranges, like their singly-charged counterparts. While, as expected, the broad 
plasmon-like structure with its maximum at about 17-18 eV is relatively insensitive 
to the charge-state of the molecule, a systematic decrease with increasing positive 
charge of the absorption cross-section between ~6 and ~12 eV is observed for each 
member of the class. 
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1 Introduction 



Polycyclic aromatic hydrocarbons [la-b] (PAHs) are a large class of conju- 
gated TT-electron systems of fundamental importance in many research areas 
of chemistry as well as in astrophysics and materials science. The carbon 
skeletons of PAHs may be considered as small pieces of graphite planes and, 
as such, they have been proposed as precursors to extended carbon networks 
such as fuUerenes and carbon nanotubes [2a-b]. PAHs are of high interest in 
environmental chemistry due to their carcinogenicity and their ubiquity as 
air pollutants produced by the combustion of organic matter [3a-b]. In the 
astrophysical context, PAHs are found in carbonaceous meteorites [4a] and 
in interplanetary dust particles [4b]. Based on the astronomical observation 
of IR, visible and UV spectroscopic features, neutral and charged PAHs are 
thought to be the most abundant molecules in space after molecular hydrogen 
and carbon monoxide [4c]. 

Oligoacenes, or simply acenes, are a subclass of catacondensed PAHs (with 
all carbon atoms on the periphery of the ring system) consisting of fused 
benzene rings joined in a linear arrangement. In their crystalline state these 
organic semiconducting materials have received particular attention in the field 
of electronics and photonics [5a-d]. Acenes and their derivatives are being 
increasingly used as active elements in a variety of opto-electronic devices 
such as organic thin-film field-effect transistors [6a-b] , hght-emitting diodes 
[7a-b], photovoltaic cells [8a-b], and liquid crystals [9]. Organic electronics 
based on functionalised acenes and heteroacenes is presently a very active 
field of research [lOa-b]. 

Since the early work by Clar [la-b] and Piatt [lla-b], there has been a 
wide-ranging interest in the electronic properties of PAHs using different spec- 
troscopic techniques, such as absorption [12a-j], electron-energy-loss [13a-d], 
fluorescence [14a-f], polarisation [15a-c], photoclcctron [16a-f], and photoion 
mass spectrometry [17a-d]. Thanks to these studies, the electronic spectra 
of neutral PAHs are known to be composed of two main regions: (i) the 
broad plasmon-like excitation peaking at ~17-18 eV, which involves n a*, 
a — > TT*, a — > <T*, and Rydberg spectral transitions, and (ii) the single-particle 
excitation part below a few cV, where the lowest energy singlet -singlet vr — tt* 
transitions occur. The four lowest transitions of neutral PAHs are usually de- 
scribed by the Clars's notation p, a, /3' [la], or ^La, ^Lb, ^Bf,, ^Ba according 
to the empirical model of Piatt [11a]. In catacondensed PAHs these transitions 
are characterised by the following intensities and oscillator strengths: p, weak, 
/ fti 0.01-0.1; a, very weak, / a; 0.001; /3, very strong, / a; 1.0; (3', medium 
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strong, / f» 0.1-1.0 (e.g., [17d]). In the case of oligoacenes, in particular, 
the transition dipole moment for the p-band hes along the short-axis of the 
molecule, while for the a and /3-bands it lies along the long-axis. In the follow- 
ing we will refer to these transitions as "short-polarised" and "long-polarised" , 
respectively. 

Charged PAHs have been the subject of extensive spectroscopic studies in 
frozen glassy organic solids [18a-b]. These experiments showed that PAH rad- 
ical cations and anions: (i) display intense optical transitions at lower energies 
than their parent molecule, and (ii) have very similar electronic spectra, in 
qualitative agreement with the particle-hole equivalence in the pairing theo- 
rem of Hiickel's theory [18a]. Currently, the interest on charged PAHs comes 
mainly from basic research in astrophysics because PAHs arc expected to ex- 
ist in space in different charge states depending on the physical conditions 
(e.g., UV flux, electron density, etc.) of the host environment [4c]. This has 
motivated a large amount of laboratory work based flrst on matrix isolation 
spectroscopy [19a-k] and, more recently, laser mass spectroscopy [20a-c] , spec- 
troscopic studies of the molecules trapped in helium droplets [21a-d], and a 
high sensivity photo-absorption technique in free jets called cavity ring-down 
spectroscopy [22a-g]. 

From an astrophysical point of view the knowledge of the electronic absorption 
spectra of PAHs in all their relevant charge states is of fundamental importance 
for our understanding of their photophysics in space. While this concerns the 
whole energetic range excitable in a typical interstellar environment, i.e., from 
the visible to the far-UV, very few experimental data are available for charged 
PAHs in this spectral range due to the limitations which are intrinsic to the 
laboratory techniques more widely used. As a part of a more extensive study 
[23a-e] towards the knowledge of the spectral properties of a large sample of 
PAHs to be modelled in astrophysical environments [24a-c] , we report in this 
paper a detailed study of the electronic absorption spectra of the five smallest 
oligoacenes naphthalene, anthracene, tetracene, pentacene, and hexacene in 
the charge states most relevant for astrophysical applications, i.e., -1, 0, +1, 
and +2. The geometries of the molecules considered are sketched in Fig. 1. 
The theoretical methods we used for both ground-state and excited-state 
calculations have been validated for the basic aromatic unit benzene in its 
neutral form, for which a large amount of spectroscopic data are available. 

There have been many systematic studies of benzenoid hydrocarbons, in- 
cluding acenes, using quantum-chemical calculations to obtain, e.g., heats 
of formation [25a], infrared spectra [25b], and C-H bond dissociation ener- 
gies [25c]. While the electronic properties of neutral acenes in the near-IR, 
visible, and near-UV spectral ranges are thoroughly characterised up to the 
sizes of hexacene [26a-k,27a-b] and some studies have been done for larger 
acenes [28a-d], a comparatively smaller amount of work exists for their cor- 
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Fig. 1. Oriented geometries of the molecules considered. Prom top to bottom: naph- 
thalene, anthracene, tetracene, pentacene, and hexacene. The molecules are sup- 
posed lying in the x-y plane, the y axis being the longer one. 

responding mono-cations [29a-h] and mono-anions [29e,29f,30a-b]. To the 
best of our knowledge, however, the electronic spectra of charged PAHs in 
the far-UV spectral domain have not been measured to date. In addition, a 
detailed study of the electronic excitation properties of PAHs in their dou- 
bly-ionised state has been missing until now. First proposed twenty years ago 
[31a-c], the possible presence of PAH dications in the interstellar medium has 
recently received further support based on the proposition that these molecules 
could contribute to the luminescence observed in the red part of the visible 
spectrum in many interstellar sources [32]. While laboratory measurements 
of the yield of optical fluorescence and phosphorescence by these species are 
needed to test this hypothesis, the determination of their electronic absorption 
spectra is motivated. 

We used density functional theory (DFT) [33a-k] and its time-dependent 
extension (TD-DFT) [34a-k], which are methods of choice for this type of 
investigations for large molecules. Since excitation energies and oscillator 
strengths within TD-DFT can be computed following two different strate- 
gies (see Sect. 2.2 for details), we used both of them to obtain: (i) the whole 
absolute photo-absorption cross-sections up to the soft x-ray region near 
30 eV, and (ii) the positions, the oscillator strengths, and the leading configu- 
rations of the lowest-lying valence tt ^ vr* short-polarised and long-polarised 
electronic transitions, falling in the visible/near-UV region, as a function 
of molecular size and charge state. In case (i), due to its numerical stabil- 
ity, we used the adiabatic local-density approximation in the parametrisa- 
tion of Perdew & Zunger [33c], which has proven to yield reliable results 
for the dynamical polarisability of conjugated molecules [23e,34g]. In case 
(ii) we used the hybrid B3LYP functional [33i] and the gradient-corrected 
BLYP functional [33d, 33e], which are widely used in the study of PAHs to 
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obtain a large number of molecular parameters, such as structures and ener- 
getics [26a,35a], vibrational spectra [35a,35b,35c], ionisation spectra [26c, 26e], 
ionisation energies [23c,26a,35f,26k], electron affinities [23b,30a,35g,35k], elec- 
tronic excitations [29d,26d,26f,26g,29f,30b,26k], the vibronic structure of ab- 
sorption spectra [26h,26i], hydrogen dissociation [29g,29h], and Jahn-Teller 
effects [35d,35e,35h,35i,35j]. 

Although the approximate exchange-correlation functionals we used have been 
developed for the electronic ground-state, they are also routinely employed in 
TD-DFT calculations, and their application usually yields accurate results for 
low-lying valence-excited states of both closed-shell [34c, 34d] and open-shell 
[34e,34f] species. It is known, however, that these functionals show the wrong 
asymptotic behaviour, decaying faster than 1/R (i.e., exponentially) for large 
distances R from the nuclei. Among the well known and documented limita- 
tions of these methods [34j,34k] are: (i) the correct description of Rydberg 
[36a-c], doubly-excited [37a-b], and charge-transfer [38a-c] excited-states, 
(ii) the failure for large, extended 7r-systems such as polyacetylene fragments 
and oligoporphyrins [39], and (iii) the systcm-sizc-dependent errors found 
for the lowest short-polarised excited states (p-bands in Clar's notation) of 
neutral oligoacenes [26g]. Furthermore, the oscillator strengths computed by 
TD-DFT are considered to be only in qualitative agreement with experi- 
ments. In a rigorous assessment of the quahty of TD-DFT molecular oscillator 
strengths, for example, they were found to be reasonable but not in quanti- 
tative agreement with reliable experimental and theoretical values for small 
molecules such as CO, N2, and CH2O [34h]. 

In spite of the above-mentioned failures of the level of theory we used, it is 
known to be sufficient to identify the most intense lowest-lying electronic exci- 
tations of neutral PAHs [26b, 26f] as well as radical ions [29d,29e,29f,30a-b,40a-c], 
which are found to match closely with the available experimental data both in 
terms of positions and intensity ratios. Despite the non-physical exponential 
fall-off of the exchange-correlation functional used, in fact, these states all 
involve excitations to and from delocalised, valence 7r-orbitals, which are not 
significantly affected by the shape of the exchange-correlation potentials in 
the asymptotic region [29f]. In particular, it has been shown that interesting 
trends exist in the vertical excitation energies and the oscillator strengths for 
homologous series of PAHs [40a,40b,40c]. For the transitions with the largest 
oscillator strength in the oligorylenes perylene, terrylene, and quaterrylene, 
for example, a net increase of the oscillator strength per unit mass of carbon 
along the series has been found [40a]. This result might have important im- 
plications in astrophysics with respect to the long-standing unsolved problem 
of the diffuse interstellar bands, about 300 unidentified absorption features 
observed in the near-UV, visible, and near-IR spectra of stars obscured by 
interstellar dust [41]. 
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Moreover, wc used the so-called delta-self-consistent-field ( ASCF) approach 
[33g], evaluating total energy differences between the self-consistent field cal- 
culations performed for the neutral and charged systems to obtain: (i) the 
vertical and adiabatic electron affinities and the first and second ionisation 
energies; (ii) the quasiparticle correction (quasiparticle energies are associated 
with the addition or removal of an electron) to the highest occupied molec- 
ular orbital (HOMO) - lowest unoccupied molecular orbital (LUMO) gap. 
This quantity is related to molecular hardness, the analogue of the band gap 
of sofids, defined as half the difference between the ionization potential and 
the electron affinity, which is a key property characterising the chemical be- 
haviour and reactivity of a molecule. The comparison between the optical gap, 
i.e., the lowest singlet-singlet excitation energy obtained with TD-DFT, and 
the quasiparticle corrected HOMO-LUMO gap enabled us to estimate the ex- 
citonic effects (due to the electron-hole interaction) in the neutral molecules. 

The paper is organised as follows. Section 2 contains the technical details for 
the calculation of the ground-state properties (Sect. 2.1) and the electronic 
absorption spectra (Sect. 2.2). The results we obtained are presented in Sect. 3 
and discussed in Sect. 4. Our concluding remarks are reported in Sect. 5. 



2 Computational details 

2.1 DFT calculations 

The calculation of the excitation energies and the electronic absorption spectra 
required the previous knowledge of the ground-state optimised geometries. For 
this part of the work we used the Gaussian-based DFT module of the nwchem 
package [42]. Geometry optimisations were performed using a basis-set with 
the smallest addition of diffuse functions, namely the 6-31-l-G* basis, a valence 
double zeta set augmented with d polarisation functions and s and p diffuse 
functions for each carbon atom. 

We used the hybrid B3LYP functional, a combination of the Becke's three pa- 
rameter exchange functional [33i] and the Lee- Yang-Parr gradient-corrected 
correlation functional [33e]. Although hybrid DFT functionals are computa- 
tionally more expensive than other exchange-correlation functionals in the lo- 
cal density or generalised gradient approximations, B3LYP results for ground-state 
properties are known to be markedly more accurate compared with experi- 
ment for a large number of systems including PAHs in general [35a, 35b] and 
oligoacenes in particular [26a,26c,26e,26k]. This is confirmed for neutral ben- 
zene, whose optimised bond-lengths we obtain at the B3LYP/6-31-I-G* level 
are rcc = 1.399 A, and rcn = 1.087 A, to be compared with the empirical 
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equilibrium (re) recommended values of, respectively, 1.3914±0.0010 A and 
1.0802±0.0020 A [43a]. Analogously, the ground-state rotational constant 
we found is ~5700 MHz, in agreement with the empirical equilibrium (Be) 
determination of 5731.73 MHz [43a]. 

Prom the structural relaxations performed for both neutral and charged sys- 
tems, we computed via total-energy differences the adiabatic electron affinities 
and the adiabatic single and double ionisation energies. At the optimised ge- 
ometry of the neutral molecule we evaluated also the vertical electron affinity 
(EAv) and the vertical first ionisation energy (lE^). This enabled us to obtain 
the quasiparticle (QP) corrected HOMO-LUMO gap of the neutral systems 
considered, which is rigorously defined within the ASCP scheme [33g] as: 

QPj^p = lEv - EAv = En+1 + En-1 - 2En, (1) 



En being the total energy of the N-electron system. We used also the following 
approximate expression: [33f|: 

QPgap ~ ~ ^JV> (2) 



where e- is the i* eigenvalue of the j-electron system. The results obtained 
using the above Eqs. (1) and (2) tend to coincide as the system gets larger and 
the orbitals more delocalised. The B3LYP/6-31-I-G* level of theory shows good 
agreement with experiments for the electron affinities of PAHs [23b,30a,35k], 
but it is known to be unable to predict their absolute ionisation energies 
with chemical accuracy (±0.1eV) [23c,26a,26k,35f]. This has been discussed 
in detail by Kadantsev et al. [26k], who concluded that a better descrip- 
tion of the electron correlation is needed to reproduce the experimental lEs. 
Switching to another computational scheme one could employ many-body 
perturbation theory in the so-called Hedin's GW approximation [44a]. This 
method, in which the QP energies are calculated from the self-energy oper- 
ator of the system (given as the product of the Green's function Q and the 
screened Coulomb interaction W), gives results in excellent agreement with 
the available experiments for many materials (sec e.g., [44b]). To assess the 
reliability of our ASCP QP-corrected HOMO-LUMO gaps, we compared them 
with the GW results obtained for the oligoacenes [28d]. The QP-corrected 
HOMO-LUMO gap of our benchmark benzene molecule is 10.59 eV (lEy = 
9.20 eV, EAy — -1.39 eV), that compares favorably with the GW results 
[28d] of 10.59 eV (first-principles calculations using Gaussian-type orbitals) 
and 10.46 eV (DFT-based tight-binding calculations), respectively, and with 
the experimental value of ~10.36 eV (lE^ = 9.24384±0.00006 [43b], EA^ = 
-1.12±0.03 eV [43c]). 
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2.2 TD-DFT calculations 



Thanks to the good compromise between accuracy and computational costs, 
compared to many-electron wavefunction-based ah initio methods, TD-DFT 
is the most widely used approach to compute the excitation energies of such 
complex molecules as PAHs [26b,26d,26f,26g,26h,26i,28b,29d,29e,29f,30b,?]. 
In this study we used two different implementations of TD-DFT in the linear 
response regime, in conjunction with different representations of the wavefunc- 
tions: 

(1) the real-time propagation scheme using a grid in real space [34g], as 
implemented in the OCTOPUS computer program [45a-d]. 

(2) the frequency-space implementation [34e] based on the linear combina- 
tion of localised orbitals, as given in the nwchem package [42]. 

In the first scheme (1) the time-dependent Kohn-Sham equations are directly 
solved in real time and the wavefunctions are represented by their discretised 
values on a uniform spatial grid. The static Kohn-Sham wavefunctions are 
perturbed by an impulsive electric field and propagated for a given finite time 
interval. In this way, all of the frequencies of the system are excited. The 
whole absolute photo-absorption cross-section <j{E) then follows from the 
dynamical polarisability «(£'), which is related to the Fourier transform of 
the time-dependent dipole moment of the molecule. The relation is: 

^{E)^^-^^{a{E)Y (3) 



where h is Planck's constant, 9{q;(£')} is the imaginary part of the dynamical 
polarisability, and c the velocity of light in vacuum. The dipole strength-function 
S{E) is related to (t{E) by the equation: 

S(E) ^ ^.W. (4, 



irie and e being respectively the mass and charge of the electron. S{E) has 
units of oscillator strength per unit energy and satisfies the Thomas-Reiche-Kuhn 

dipole sum-rule = J dE S{E), where is the total number of electrons. 

The great advantage of obtaining the whole response at once is particularly 
useful for astrophysical applications, for which the whole absorption spectrum 
is needed. 

In the most widely used frequency-space TD-DFT implementation (2), based 
on the linear response of the density-matrix, the poles of the linear response 
function correspond to vertical excitation energies and the pole strengths to 
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the corresponding oscillator strengths [34b] . With this method computational 
costs scale steeply with the number of required transitions and electronic ex- 
citations are thus usually limited to the low-energy part of the spectrum. 
Prom a computational point of view the advantages of the real-time propa- 
gation method are discussed, e.g., in Ref. [45d]. On the other hand, the main 
drawbacks of the real-time approach are that: (i) no information is given on 
dipole-forbidden singlet-singlet and singlet-triplet transitions, and (ii) one 
does not obtain independent information for each excited state, such as its 
irreducible representation of the point group of the given molecular system, 
and the description of the excitations in terms of promotion of electrons in an 
orbital picture. 

We performed the OCTOPUS calculations in the local-density approximation, 
with the exchange-correlation energy density of the homogeneous electron gas 
[46a] parametrised by Perdew & Zunger [33c] . The ionic potentials are replaced 
by norm-conserving pseudo-potentials [46b] . We used a grid spacing of 0.3 A 
and determined the box size by requiring each atom to be at least 4 A away 
from its edges. We furthermore added a 1 A thick absorbing boundary, which 
quenches spurious resonances due to standing waves in the finite simulation 
box used to confine the molecules [34g,45c] . We used a time integration length 
T=20 h/eV, corresponding to an energy resolution of ^/T=0.05 eV. For the 
numerical integration of the time evolution we used a time step of 0.002 h/eV, 
which ensured energy conservation with good numerical accuracy. 

The TD-DPT calculations with nwchem were performed at the same level 
B3LYP/6-31-I-G* used to obtain the ground-state geometries. Although basis 
set convergence is not yet expected at the level we used, our results for the 
neutral systems are almost coincident with the ones obtained in Ref. [26k] 
using the larger 6-311-I--I-G** basis, which is supplemented with a third layer 
of valence functions and includes polarisation and diffuse functions on both 
carbon and hydrogen atoms. We thus believe our theoretical predictions to 
be sufficiently accurate for the purposes of this work. In the case of neutral 
benzene, we predict the strong vr — tt* transition ^Aig — s>^Eiu at 6.96 eV with 
an oscillator strength of 1.22, in good agreement with the measured band 
position in vapour-phase of 6.94 eV with an /-value of 1.2 [43d]. 

In order to assess the choice of this specific exchange-correlation functional, at 
the B3LYP optimised geometries, we used also the gradient-corrected BLYP 
functional [33d, 33e], the same approach used in previous studies of PAH 
ions [29d,40a,29f,30b]. With both methods we restricted ourselves to the 
first 20 singlet-singlet roots. Since we are interested in the behaviour of the 
lowest-lying permitted in-plane long-polarised and short-polarised electronic 
transitions as a function of charge state and molecular size, in the following 
we report only the first five electronic transitions. The complete set of elec- 
tronic excitation energies and oscillator strengths computed at both B3LYP 
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and BLYP levels, including also the optically inactive ones, are available in 
our online database of the computed spectral properties of PAHs [23d] . 



3 Results 

3.1 Static properties 

Geometry optimisations with nwchem were performed using tight conver- 
gence criteria, that are specified by maximum and root mean square gradient 
thresholds of 1.5-10"^ and 1.0-10^^ atomic units, respectively, and maximum 
and root mean square thresholds of the Cartesian step respectively of 6.0-10"^ 
and 4.0-10"^ atomic units. According to previous studies [29b, 29d] the lifting 
of the molecular symmetry D2h of the neutral molecules is not expected to lead 
to optimised geometries with lower symmetry for the corresponding charged 
species. We indeed confirmed that the structural parameters obtained for 
naphthalene, anthracene, and tetracene in the charge-states and +2, ig- 

noring any apparent symmetry and adopting the D2h symmetry, are coincident 
within numerical errors. We therefore assumed the above D2h constraint for all 
subsequent calculations in the paper. Our B3LYP/6-31-I-G* geometry optimi- 
sations for the neutral molecules give structural parameters in good agreement 
with those previously published [26a,26e,26k]. In particular, although a larger 
basis was used in Ref. [26k], the two sets of results are almost coincident and 
compare fairly well with the available x-ray data. We do not discuss here the 
changes of the single bond lengths and bond angles occurring in the charged 
species, compared to the corresponding neutral ones. Depictions of the struc- 
tures of each molecule considered, in which internal coordinates are shown and 
compared, are given as supplementary material to the paper. Instead, Fig. 2 
presents in a collective way the structural variations relative to the neutral 
molecules, expressed in terms of the percentage variations of the rotational 
constants A and B. These latter quantities are proportional to the inverse 
of the principal momenta of inertia corresponding to the in-planc short and 
long-axis of the molecules, respectively: A=(h/87r^c)Ishort) B=(h/87r'^c)Iiong- 
Note that in Fig. 2 we omitted the results corresponding to naphthalene an- 
ion, since gas-phase naphthalene is unable to bind an additional electron 
in its LUMO state [43c]. The ground-state optimised geometries of all of 
the molecules considered, both Cartesian and internal coordinates, are freely 
available in our online database of the computed spectral properties of PAHs 
[23d]. 

All neutral and singly-charged species were computed as singlet and doublet, 
respectively, while for dications we computed both their singlet and triplet 
ground-states. The adiabatic and vertical values of electron affinities and sin- 
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Fig. 2. Percentage variation of the rotational constants A (top panel) and B (bottom 

panel) relative to the neutral counterparts for anions (asterisks), cations (triangles), 
and dications (diamonds), as a function of molecular size. Wc omit the entries for 
naphthalene anion, which is known not to form a stable anion in the gas-phase. 

gle and double ionisation energies as obtained via total energy differences are 
given in Table 1 and compared with the available experimental data [47,31c]. 
As shown in Table 1, for all of the five molecules considered in their dou- 
bly-ionised state, our calculations predict the total energy of the singlet state 
to be lower than that of the triplet state by ~0. 5-1.0 cV. Figure 3 displays 
the computed electron affinities and first and double ionisation energies as a 
function of size, and compares them with the available laboratory data. 
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Tabic 1 

Adiabatic and vertical values (in parentheses), all data in eV, of electron affinities and single and double ionisation energies of the 
oligoacenes considered in this work (C4n+2H2n+4, n=2, 3, 4, 5, 6) as obtained through total energy differences at the B3LYP/6-31+G* 
level. For comparison we list also the experimental adiabatic electron affinities and adiabatic single ionisation energies taken from the 
NIST Chemistry WebBook [47], as well as the adiabatic second ionisation energies photon-impact measurements [31c]. 



Number 


Electron affinity 


First ionisation energy 


Dication 


Double ionisation energy 


of cycles 


Ad. (Vert.) 


Exp. 


Ad.(Vert.) 


Exp. 


state 


Ad.(Vert.) 


Exp. 


2 


-0.26(-0.38) 


-0.20±0.05" 


7.80(7.89) 


8.144±0.001^ 


singlet 
triplet 


20.99(21.35) 
21.45(21.57) 


21.5±0.2 


3 


0.53(0.43) 


0.530±0.005^ 


7.02(7.09) 


7.439± 0.006'^ 


singlet 
triplet 


18.70(18.95) 
19.70(19.80) 




4 


1.08(1.00) 


1.067±0.043^ 


6.49(6.55) 


6.97±0.05-^ 


singlet 
triplet 


17.15(17.34) 
17.96(18.11) 


18.6± 0.2 


5 


1.48(1.41) 


1.392±0.043^ 


6.12(6.16) 


6.589±0.00F 


singlet 
triplet 


16.03(16.18) 
16.67(16.80) 


17.4±0.2 


6 


1.78(1.72) 




5.83(5.87) 


6.36±0.02-^ 


singlet 
triplet 


15.18(15.30) 
15.68(15.78) 





" Extrapolated from the EAs of naphthalene-water clusters determined via photoelectron spectroscopy [16e]. 

From laser threshold photoelectron spectroscopy [16c]. 
^ From photodetachment photoelectron spectroscopy [16d]. 

From two-laser photoionisation supersonic jet mass spectrometry [17a]. 
^ Estimated from gas-phase electron attachment free energies with the electron-transfer equilibria technique [48]. 
^ From gas-phase photoelectron spectroscopy [16b]. 
^ From high-resolution gas-phase photoelectron spectroscopy [16f]. 
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Fig. 3. Computed adiabatic ionisation energies and electron affinities of tfie stud- 
ied oligoacenes as a function of size. The corresponding experimental values are 
represented by the filled symbols (see Table 1). 

3.2 Excitation energies and electronic absorption spectra 



Since its first applications, the real-time TD-DFT method in real space was 
proven to give good results for neutral benzene [34g,45c], compared with the 
experimental spectrum recorded in the energy range 6-35 eV [12e]. Applying 
this same approach to a large sample of PAHs, we already showed [23a, 24b] 
our results to be in good agreement up to photon energies of about 30 eV 
with the experimental data obtained for a few neutral PAHs with the syn- 
chrotron radiation facility of SUPERACO [12h,12i]. The comparison between 
the theoretical spectra obtained with the OCTOPUS code and the experimental 
photo-absorption cross-sections of neutral anthracene (C14H10) and benzene 
in the gas-phase are shown in Figs. 4 and 5. The latter experimental spectra 
are in good agreement with the ones obtained with the synchrotron radia- 
tion from the electron accelerator DESY for anthracene [12g], and benzene 
[12e], respectively. Since the theoretical spectrum is averaged over the three 
X, y, and z polarisations, each single contribution is also shown in Figs. 4 and 
5. As examples, the spectra computed for tetracene (C18H12) and pentacene 
(C22H14) in the four charge states considered are given in Figs. 6 and 7. The 
spectra for all of the other molecules under study can be found in our online 
database [23d]. Figure 8 shows the comparison between the integrated values 
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Fig. 4. Comparison between the computed (solid black line) photo-absorption 
cross-section a{E) of neutral anthracene (C14H10) and the corresponding gas-phase 
absorption spectrum (dotted line, taken from Refs. [12h,12i]). The contributions 
corresponding to polarisations along the x-axis (in-plane short), y-axis (in-plane 
long), and z-axis (out-of-plane) , are marked in gray, dark gray, and light gray, 
respectively. Units are megabarns, 1Mb = 10"-*^^ cm^. 

in the range 6-12 eV of the individual dipole strength-functions S{E) (see 
Eq. 4) divided by the number of carbon atoms in each molecule as a function 
of molecular size. 



The first few permitted electronic transitions of each molecule, as obtained at 
the B3LYP/6-31+G* and BLYP//B3LYP/6-31+G* levels with the TD-DFT 
frequency-space implementation of nwchem are reported in Tables 2-6, and 
compared with the available experimental data that we could find in the litera- 
ture. We deliberately omitted the large amount of photoelectron data available 
for neutral species [16a-f|. The use of such data for spectral assignments of the 
so-called Koopmans transitions of radical cations has been already discussed 
many times [29c,29d,29e,29f]. Electronic excited states are classified under 
the point-group and the ground-state symmetry is specified for each 
charge-state. As to the character of the excited electronic states, we analyse 
the nature of the corresponding transitions in terms of the occupied and virtual 
molecular orbitals that have been interchanged between the ground and the 
excited electronic states [34j] . The above description is computationally conve- 
nient and straightforward for states which are well described with only one or 
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Fig. 5. Same as Fig. 4 for neutral benzene. Note that, due to the symmetry of the 
molecule, the curves corresponding to x and y polarisations are obviously coincident. 

two significantly contributing "excited" Slater determinants [34j]. This is the 
case for the electronic excited states reported in Tables 2-6. The use of more 
sophisticated and physically more appealing ways to obtain information about 
an electronic transition, such as difference density or attachment / detachment 
density plots [34j] , is outside the scope of the present work. We used the same 
notation as in Refs. [29d,29f] where the tt orbitals are numbered in the order 
of increasing energies and 7r_i, 7ro(7rQ), and tt* denote the highest doubly oc- 
cupied TT orbital, the singly occupied (unoccupied) tt orbital, and the lowest 
doubly unoccupied tt orbital, respectively. 

Figures 9 and 10 display, respectively, the B3LYP/6-31-I-G* positions and 
the corresponding oscillator strengths of the lowest in-plane short-polarised 
(p-bands in the neutral molecules) and long-polarised (a-bands in the neutral 
molecules) electronic transitions as a function of the number of benzene units 
and the charge state of the molecule. 



3.3 Quasiparticle-corrected HOMO-LUMO gap of the neutral species 

For each of the neutral molecules considered at the B3LYP/6-31+G* level, 
we report in Table 7 the comparison between the HOMO-LUMO gap Eg^^p 
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Table 2 

Singlet-singlet lowest-lying permitted excitation energies (in eV) of naphthalene 
anion, neutral, cation, and dication as obtained via frequency-space TD-DFT. The 
corresponding oscillator strengths are given in parentheses. Polarization of 
the bands are denoted according to Fig. 1 as x, y, and z for in-plane short, in-plane 
long, and out-of-plane polarized bands, respectively. The description of each 
excitation is based on the B3LYP results and is given in terms of the 
occupied and virtual molecular orbitals contributing significantly to it 
[29d,29f|. We report also the available experimental data for comparison. 



State (pol.) 


Excitation 


BLYP 


B3LYP 


Exp. 






Anion (ground-state ^Bsg) 




1 'B2„ (z) 


TTo 7r| 


0.66(0.001) 


0.99(0.001) 




1 'Bi„ (y) 


TTO TTy 


1.48(0.022) 


1.58(0.034) 


~1.5" 


2 ^B2u {z) 


TTo TT^ 


1.85(0.005) 


2.15(0.006) 




2 'Bi, (y) 


TTo (72 


2.30(0.071) 


2.49(0.069) 


~2.5" 


1 2A„ (x) 


7r_l TT^ 


2.73(0.001) 


2.89(<0.001) 





1 ^B3„ {x,p) 
1 ^B2„ {y,a) 
2^B2„ (y,/3) 

2 iB3„ 

l^Bi 



7r-2 
7r-2 



X) 



Neutral (ground-state -"^A^) 
7r_i ^ ttJ 4.06(0.045) 4.36(0.061) 

4.20(<0.001) 4.44(<0.001) 
5.62(1.155) 5.85(1.260) 
7r_2 ^ TT^ 5.80(0.146) 6.08(0.199) 

7r-2 ^ TT^ 5.77(0.008) 6.24(0.016) 



TT*, 7r_i 
TT*, 7r_i 



^2* 



TTo 



4.45(0.102)^ 4.45(0.109)^ 4.44^^ 
3.97(0.002)^ 4.0^ 3.98'^ 
5.89(1.3)^ 5.89(1.3)^ 5.86'' 
6.14(0.3)^ 6.0^ 6.09'^ 



l'B25 

l^Bs, 

2 'B33 
2'B2, 

3 ^623 



y) 

2/) 



7r-2 
7r-3 

TTo - 
TTo - 



Cation (ground-state ^A^) 
2.15(0.042) 2.14(0.053) 



TT, 



2.77(0.006) 
3.51(0.051) 
3.74(0.012) 
4.30(0.008) 



2.98(0.006) 
3.59(0.064) 
3.91(0.032) 
4.61(0.006) 



1.84<^, 1.85(0.052)^ 1.84(0.011)-^ 
2.72'^, 2.72(0.010)^ 2.69(0.001)-^ 
3.29'^, 3.25(0.016)-^ 
4.02*^, 4.03(0.024)^ 
4.49'=', 4.52(0.060)-^ 



1 ''^B2u 

l^Bi 



2 ^B2u 



1 ^B, 



2 ^Bi 



y) 

x) 

y) 

z) 



7r-2 
7r_4 
7r-3 

TT-l 



4 



Dication (ground-state ^A^) 
2.85(0.061) 3.00(0.090) 
2.71(<0.001) 3.25(<0.001) 
3.32(0.062) 3.67(0.090) 
5.06(0.833) 5.28(0.865) 
4.95(0.001) 5.58(0.002) 



"Absorption in glassy organic solid [18a]. ^Gas-phase absorption [12c]. 
'^Electron-energy-loss spectroscopy [13b]. Absorption in neon-matrix [19d]. 
^Gas-phase absorption [20b]. -^Absorption in argon-matrix [19b]. 
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Energy (eV) 



Fig. 6. Computed photo-absorption cross-section a{E) of tetracene (C18H12) in neu- 
tral (black), anionic (dark gray), cationic (gray) and dicationic (light gray) charge 
state, as obtained with the real-time real-space implementation of TD-DFT. 

obtained as difference of Kohn-Sham eigenvalues, the excitation energy of the 
HOMO-LUMO transition E™-^^'^ as given by TD-DFT, and the corre- 
sponding experimental value Eg^P (p-bands in Tables 2-6). In the same table 
we compare the results of Eqs. (1) and (2) with the DFT-based tight-binding 
GW data of Ref. [28d], QPg^^p^"^^, as well as with the corresponding exper- 
imental value obtained as the difference between the experimental EAs and 
first lEs given in Tabic 1, QPgap = lEcxp-EAcxp- The theoretical exciton bind- 
ing energy Ebind is estimated through the difference QPgg^p-E™~°^'^, which is 
compared with its corresponding experimental value QP|ap - E|ap- All of these 
quantities are displayed in Fig. 11 as a function of molecular size. 
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Fig. 7. Same as Fig. 6 for pentacene (C22H14). 



0.5 














0.4 


- )K 


yd 




































; + 












0.3 




+ 


+ 
















+ 


+ 






': A 


A 


A 


A 


A 




0.2 


\ 
















0.1 














2 


3 


4 


5 


6 








Number 


of 


cycles 







Fig. 8. Comparison between the integrated values in the range 6-12 eV of the indi- 
vidual dipole strength-functions S (E) (see Eq. 4) divided by the number of carbon 
atoms for the oligoacenes anions (asterisks), neutrals (crosses), cations (triangles), 
and dications (diamonds) considered, as a function of molecular size. 
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Fig. 9. Calculated positions (B3LYP/6-31+G*) of the lowest short-polarised 
(crosses, p-bands in the neutral species) and lowest long-polarised (diamonds, 
a-bands in the neutrals) electronic transitions in the five oligoacenes considered 
as a function of the size and charge-state of the molecule. 
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Fig. 10. Same as Fig. 9 for the corresponding oscillator strengths. 
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Fig. 11. Top panel: optical gap E^^,^-"^ ^ as obtained via TD-DFT at the 

B3LYP/6-31+G* level (open triangles) and ASCF QP-corrected HOMO-LUMO 
gap QPgap computed via Eq. (1) (open diamonds), as a function of molecular size. 
Bottom panel: exciton binding energy Ebind estimated as QPgap - E™"^^""" (open 
squares). In both figures, the corresponding experimental values are represented by 
the filled symbols (see Table 7). 
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Table 3 

Same as Table 2 for anthracene in its -1,0, +1, and +2 charge states. 



State (pol.) 


Excitation 


BLYP B3LYP 


Exp. 






Anion (ground-state ^Bi„) 




1 (z) 


TTo TT^ 


1.23(0.001) 1.56(0.001) 


— 


1 (y) 




1.83(0.110) 1.90(0.137) 


~1.7« 


2 (z) 


TTO ^ TT^ 


1.88(0.001) 2.17(0.002) 


— 


1 ^B2g (x) 


TT-l ^ TT^ 


2.12(0.008) 2.20(0.012) 


~2.4" 


2 ^B2g (X) 


4- 

TTo TT Ji 


2.78(0.008) 3.01(0.010) 








Neutral (ground-state -"^A^) 




1 ^Bsu {x, p) 


TT-l ^ TTi 


2.92(0.039) 3.21(0.058) 


3.27(0.1)^ 3.45^ 3.43^* 


1 'B2„ {y,a) 


7r_2 ^ TT^, 7r_i - 


^ ttI 3.60(<0.001) 3.85(<0.001) 


3.47-3.60-3.45% 3.84-^ 


in 1 / \ 

2 ^B3„ (x) 


TT-l ^ TTg 


4.59(<0.001) 5.06(<0.001) 




^ ■D2m [Vj P) 


-77- _ V TT^ TT I 

71-2 ^ T^i , 7r_l — 




A S4('9 98"!^ 94"= 




7r_l > TTy 










Cation (ground-state ^B2g) 




1 (y) 


7r-2 ^ < 


1.87(0.085) 1.93(0.108) 


1.71", 1.73(0.076)5', 1.75'* 


1 ^Biu (x) 


TTo ^ TT^ 


2.21(0.004) 2.32(0.007) 


2.02(0.018)3' 


2 ^Bi„ (x) 


7r_4 7r5 


2.91(0.044) 3.17(0.058) 


2.83 , 2.90(0.026)3 


2 (y) 


7r_i ^ ttJ 


3.22(0.017) 3.39(0.054) 


3.52(0.104)3 


3 (y) 


TTo ^ TT^ 


3.69(0.016) 4.02(0.010) 


3.95(0.187)3 






Dication (ground-state ^Ag) 




1 'B2„ (y) 


7r_2 — tt\ 


2.31(0.139) 2.45(0.206) 




1 'B3„ (rr) 


7r_3 ^ ttJ 


3.02(0.052) 3.39(0.078) 




2 ^B2„ (y) 


TT-l 1^2 


4.38(1.368) 4.63(1.444) 




1 'Bi„ (z) 


7r_8 ^ ttJ 


3.92(<0.001) 4.66(<0.001) 




2 'B3„ (x) 


7r_4 — ^ 7r2 


4.96(0.002) 5.59(0.001) 





"Absorption in glassy organic solid [18a]. ^Absorption in n-heptane solution [lib]. 
■^Gas-phase absorption [12g]. ''Fluorescence jet spectroscopy [14b]. "^Magnetic circu- 
lar dichroism measurements in solvents [15a,15b,15c]. -^Two-photon absorption in 
solution [12d]. 3 Absorption in argon-matrix [19e]. '*Gas-phase absorption [22d]. 
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Table 4 



Same as Talilo 


2 for (ct raccuo in it.: 


-1.0. + ;] , and +2 ( 


■liarg(^ siat(^s. 




State (pol.) 


Excitation 


BLYP 


B3LYP 


Exp. 






Anion (ground-state ^B^g) 




1 (x) 


TT-l ^ TT^ 


1.53(0.009) 


1.60(0.014) 


1.69" 


1 'Bi„ (y) 


TTo ^ TtJ 


1.68(0.164) 


1.77(0.209) 


1.50« 


1 'B2„ (z) 


TTo TT^ 


1.68(0.001) 


2.06(0.001) 


1.9P 


2 2Bi„ (y) 


TTo ^ TT^o 


2.40(<0.01) 


2.62(<0.001) 




2 ^B2u (z) 


TTO TT^i 


2.58(0.001) 


2.91(0.001) 











Neutral (ground-state ^Ag) 




1 ^B3„ {x,p) 


7r_i - 




2.16(0.031) 2.45(0.049) 


2.62(0.08)^ 2.60(0.11)^ 2.72'^ 


1 ^B2„ (y,a) 


7r_3 ^ ttJ, 


7r_i -> 


ttJ 3.21(0.001) 3.47(0.002) 


3.12<^ 


2 ^B3„ {x) 


TT-l - 




3.98(<0.001) 4.59(<0.001) 




2^B2„ (y,/3) 


7r_3 ^ TT^, 


7r_i -> 


TT^ 4.34(2.317) 4.62(2.691) 


4.55(1.85)^ 4.50(1.75)^ 


3 'B3„ {x) 


7r_4 - 




4.35(0.010) 4.92(<0.001) 










Cation (ground-state ^A^) 




1 'B2, (y) 


TT-l - 




1.59(0.129) 1.70(0.169) 


1.43", 1.43^ 


1 'B3g (x) 


TTo ^ TtJ 


1.58(0.008) 1.70(0.012) 


1.65«, 1.66^ 


2 2B3g (x) 


TT— 4 




2.67(0.029) 3.03(0.041) 


3.14«, 3.16^ 


2 'B2, (y) 


7r-2 




2.81(0.025) 3.05(0.084) 




3 'B2, (y) 


TT — 5 




3.17(0.005) 3.52(0.005) 










Dication (ground-state ^A^) 




1 'B2„ (y) 


TT-l - 




1.97(0.226) 2.11(0.335) 




1 'B3„ (x) 


7r_4 - 




2.80(0.042) 3.16(0.065) 




2 ^B2„ (y) 


7r-5 - 




3.65(0.093) 4.05(0.626) 




1 'Bi„ (z) 


7r_8 - 




3.46(<0.001) 4.11(<0.001) 




3 ^B2„ (y) 


7r_2 - 




3.86(1.709) 4.16(1.353) 





"Absorption in glassy organic solid [18a]. ''Absorption in n-heptane solution [lib]. 
Absorption in benzene solution [12a]. ''Absorption in gas-phase [14a]. '^Absorption 
in argon-matrix [19c, 19g]. 
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Table 5 



Same as Table 2 for pcntacene in its -1,0, +1, and +2 charge states. 


State (pol.) 


Excitation 


BLYP 


B3LYP 


Exp. 








Anion (ground-state 


'Bi„) 




1 'B2„ (x) 


7r_i - 


TTo 


1.10(0.008) 


1.16(0.013) 


1.06", 1.37^^ 


1 'B3„ (y) 


''0 


> TT* 
^ ^1 


1.50(0.224) 


1.60(0.286) 


1.40°', 1.42^= 


1 ^A, {z) 




^ "3 


1.96(0.001) 


2.35(0.001) 




2 ^A, {z) 


TTn — 


^ "8 


2.53(<0.001) 


2.81(0.001) 


— 


2 'Bag (y) 


TT— 1 - 


— »■ /I4 


2.65(0.016) 


2.86(0.118) 


2.82°^ 








Neutral (ground-state ^Ag) 




1 ^B'^,, fx, z*) 


TT-l - 




1.63(0.023) 


1.91(0.041) 


2.12(0.08)^ 2.28*=, 2.31*^ 


1 "'^Bo,, fy, a) 


7r_3 — ^ TT*, 


7r_i — 


7r| 2.94(0.004) 


3.21(0.005) 


3.73^= 


2 ^Bsn (x) 


TT /I - 


"1 


3.33(<0.001) 


3.95(<0.001) 




2 iB2„ (y,/3) 


" — O ' " 1 ) 


TT 1 — 


^ vr^ 3.93(2.706) 


4.24(3.346) 


4.00(2.20)^ 4.40^= 


3 ^B3„ (x) 


TT 1 - 


/I4 


3.66(0.021) 


4.26(0.003) 


— 








Cation (ground-state 


'B2g) 




1 ^Bi„ fx) 


TTo - 




1.17(0.007) 


1.25(0.012) 


1.27^^ 


1 (y) 


TT-l - 




1.42(0.182) 


1.50(0.238) 


1.31"^ 


2 (y) 


7r_2 - 




2.61(0.034) 


2.80(0.122) 


2.92^= 


2 2Bi„ (x) 


7r_4 - 




2.60(0.019) 


2.85(0.013) 




2 'Bi„ (x) 


7r_4 - 




2.85(0.001) 


3.10(0.022) 










Dication (ground-state ^A^) 




1 'B2„ (y) 


TT-l - 




1.73(0.324) 


1.86(0.483) 




1 'B3„ (x) 


7r_4 - 




2.61(0.029) 


2.98(0.050) 




2 ^Ba. (y) 


7r-5 - 




3.17(0.080) 


3.57(0.240) 




2 ^Bsn (x) 


7r-3 - 




3.20(0.011) 


3.66(0.001) 




3 ^B2u (y) 


7r-2 - 




3.44(2.022) 


3.75(2.186) 





"Absorption in glassy organic solid [18a]. ''Absorption in n-heptane solution [lib]. 
■^Absorption in neon-matrix [29e]. ''Absorption in gas-phase [12j]. 
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Tabic 6 



Same as Table 2 for hcxaccnc in its -1,0, +1, and +2 charge states. 



State (pol.) 


Excitation 


BLYP 


B3LYP Exp. 




Anion (ground-state ^Ba^) 




1 {x) 




0.78(0.006) 


0.83(0.011) — 


1 'Bi, (y) 


TTo -K* 


1.34(0.284) 


1.43(0.370) — 


2 {x) 


7r_2 ^ TT^ 


2.41(<0.001) 


2.47(<0.001) — 


1 'B2„ {z) 


TTo TT^ 


2.14(<0.001) 


2.56(0.001) — 


2 2Bi„ (y) 


7r_l ^ TT^ 


2.51(0.043) 


2.71(0.149) — 



Neutral (ground-state ^Ag) 

1 iB3„ (.T,p) 7r_i^7r^ 1.24(0.017) 1.51(0.034) 1.90' 

liB2„(y,a) 7r_3 ^ TT^, 7r_i ^ ttJ 2.75(0.009) 3.02(0.010) 2.80' 

2 ^Bs^ (x) 7r_4 ^ 7rt 2.78(0.001) 3.39(<0.001) — 

3 ^B3„ (x) 7r_i ^ 7r| 3.07(0.028) 3.69(0.010) — 



4 ^B3„ (x) 7r-2 ^ 7r| 3.35(0.017) 3.82(0.061) 

Cation (ground-state ^A„) 



1 'B3, 


[x) 


TTo - 




0.84(0.006) 


0.90(0.010) 


1^B2, 


(y) 


TT-l 




1.27(0.246) 


1.34(0.317) 


2 ^B3g 


(x) 


7r_3 




2.30(<0.001) 


2.40(<0.001) 


2 ^B2g 


(y) 


7r_2 




2.43(0.046) 


2.63(0.167) 


3 ^B3g 


(x) 


7r_4 




2.52(0.007) 


2.92(0.021) 



Dication (ground-state ^A^) 



1 'B2„ 


(y) 


TT-l 




1.55(0.440) 


1.67(0.652) 




1 B3„ 


(x) 


7r_4 




2.43(0.009) 


2.83(0.031) 




2 ^B3.u 


(.t) 


7r-3 




2.66(0.018) 


3.05(0.001) 




2 "'^B2u 


(y) 


7r-5 




2.82(0.062) 


3.21(0.130) 




3 ^B2« 


(y) 


7r-2 




3.10(2.203) 


3.44(2.642) 





"Extrapolated from solution spectra to the gas-phase (see compilation in Ref. [26k]). 
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Tabic 7 

Comparison between the B3LYP/6-3H-G* results (all values in eV) for the HOMO-LUMO energy gap E^^^p obtained as difference of 
Kohn-Sham eigenvalues, the excitation energy of the HOMO-LUMO transition as given by TD-DFT, E™-°^'^, and its corresponding 
experimental value Egap (jo-bands in Tables 2-6). The QP-gap evaluated using Eqs. (1) and (2) is compared with the DFT-based 
tight binding GW data of Rcf. [28d], QPg^^p"^"*^^ , as well as with the experimental value QPgap obtained as the difference between 
the experimental EAs and first lEs given in Table 1. The theoretical and experimental exciton binding energy Ebind are estimated as 
QPgap"-^gap~°^'^ QPgap - Egap, respectively. For comparison, the results of Ref. [26k] are reported within parentheses. 



n 


gap 


pTD-DFT 

gap 


TTiCXp 

^gap 


QPgap 


QPgap 


npDFT-GW 
^-■^gap 


p,pcxp 
SjJ^gap 


Ebind 


■pexp 
-^bind 


2 


4.74(4.75) 


4.36(4.35) 


4.45 


8.27(8.29) 


8.12 


8.0 


8.34 


3.91 


3.89 


3 


3.54(3.55) 


3.21(3.21) 


3.45 


6.66(6.60) 


6.58 


6.6 


6.91 


3.45 


3.46 


4 


2.74(2.75) 


2.45(2.44) 


2.72 


5.55(5.56) 


5.50 


5.5 


5.90 


3.10 


3.18 


5 


2.19(2.19) 


1.91(1.90) 


2.31 


4.75(4.76) 


4.72 


4.8 


5.20 


2.84 


2.89 


6 


1.78(1.78) 


1.51(1.50) 


1.90 


4.15(4.16) 


4.13 


4.3 




2.64 





4 Discussion 



Figure 2 shows that structural variations between charged species and their 
respective neutral counterparts display the same well-defined trend as a func- 
tion of molecular size: a general decrease for Aid (with the possible exception 
of anthracene anion), and an increase of Brci is observed for all charge-states 
considered. More specifically, the largest structural changes are observed for 
the anions, a consequence of the strongly antibonding character of the LUMO 
of the neutral counterpart. The variations in both cations and dications are 
nearly equal and remain small. A shortening along the short-axis (Arei > 0) 
compared to the neutral molecules is always observed. The anions appear 
to be primarily distorted along the long axis showing a lengthening along 
it (Brei < 0) for all sizes. These findings agree with a Hartree-Fock study of 
naphthalene and anthracene cations [25b], and a DFT study of anthracene an- 
ion performed with both hybrid and gradient-corrected exchange-correlation 
functional [30a]. 

Table 1 and Figure 3 confirm for the oligoacenes the good agreement found 
for the whole class of PAHs between the B3LYP/6-31-I-G* electron affinities 
and the available experimental data. The differences between computed and 
observed first ionisation energies are of the same order of magnitude as in pre- 
vious analyses [26a,26k], and seem to increase sligthly at increasing molecular 
size (from ~4% for naphthalene to ~10% for hexacene). Double ionisation en- 
ergies are found to change more rapidly along the series than single ionisation 
energies, and the relative errors of our calculations, in comparison with the 
only three experimental data available, increase from ~2% for naphthalene to 
~8% for anthracene and tetracene. 

As shown in Figs. 4 and 5, the real-time real-space TD-DFT method pro- 
vides results in very good agreement with the available experimental data 
for neutral species up to about 30 eV. The broad plasmon-like excitation 
peaking at about 17.5 eV, which involves tt a*, a vr*, a — > o"*, and 
Rydberg spectral transitions, is well reproduced both in position and width. 
However, as previously discussed in the case of neutral benzene [34g] , the use 
of the finite simulation box and the absorbing boundary at its edges does not 
give a satisfactory treatment of continuum efi^ects producing spurious struc- 
tures. The contribution of the three possible polarizations to the total absorp- 
tion cross-section varies considerably. While the low-energy part is due to 
in-plane-polarised electronic transitions {x and y axes), the contribution of 
the 2;-axis perpendicular to the plane of the molecule is significant only above 
a few eV. These features are found to be common to all PAHs [23a-e]. In the 
case of oligoacenes, due to their special symmetry, the strongest absorption {j3 
band in the neutral molecules) corresponds to long-axis polarization, which 
can be simply understood as the classical resonance in a conducting rod [11a]. 
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Prom Figs. 6 and 7 it is seen that the broad plasmon-hke stnictTirc with 
its maximum at 17-18 eV is relatively insensitive to the charge-state of the 
molecule. On the other hand, the charge-state of the molecule shows up in the 
low-energy range and, interestingly, the onset of this broad absorption moves 
blue-ward and becomes steeper with increasing positive charge. As shown 
in Fig. 8 this translates into a systematic decrease with increasing positive 
charge of the absorption cross-section in the energy gap between the tt and a 
plasmon-like structures (between ~6 and ~12 eV). This same behaviour has 
been shown for a larger sample of PAHs [23c] . In addition, as expected, the 
observed scatter between the different charge-states decreases with increasing 
molecular size. 

Concerning the visible-UV part of the spectrum, as reported in Tables 2-6, we 
checked both the sensitivity of our TD-DFT calculations to the use of different 
exchange-correlation functionals, and their reliability in comparison with the 
available experimental data. General agreement is found between BLYP and 
B3LYP which yield the same ordering of the strongest dipole-allowed excited 
states. On the average, the BLYP energies are found to be systematically 
smaller by 0.1-0.6 eV compared to the corresponding B3LYP results Our data 
confirm the spectral assignment done in previous studies for neutral [26k], and 
singly-charged species [29d,29f]. From the more accurate experimental data 
reported for the anions, neutrals, and cations we found that the mean relative 
deviation of the B3LYP functional is of the order of 6%, compared to the 
mean relative deviation of about 7% given by BLYP. We therefore consider 
only the former set of results in the following analysis. 

The frequency-space implementation of TD-DFT enabled us to gain some 
insight into the nature of the first few electronic excitations. Focusing only 
on dications, reported here for the first time, we find that the lowest in-plane 
long {y) and short-axis {x) polarised bands correspond, respectively, to the 
HOMO-l^LUMO (7r_2 ^ tt*) and H0M0-2^LUM0 (7r_3 ^ transi- 
tions for naphthalene and anthracene, and HOMO^LUMO (/T_i ti\) and 
H0M0-3^LUM0 (7r_4 tt^) for tetracene, pentacene, and hexacene. Fig- 
ures 9 and 10 show interesting trends as to the behaviour of the lowest-lying 
electronic transitions as a function of the size of the molecule and its charge 
state. The positions of both in-plane short and long-polarised excitations 
are found to shift to lower energies with molecular size for all charge states 
considered (Fig. 9). We observe the well-known similarities in the electronic 
absorption spectra between anionic and cationic PAHs, as well as the system- 
atic shifts in band position when going from the cation to the anion [18a,29f]. 
The sign and the magnitude of these shifts, attributed to the different effect 
of the (T-electrons in both ions [18a], are reproduced by TD-DFT for the most 
intense bands [29f], i.e., the lowest-lying y bands in oligoacenes. For example, 
the measured blue-shifts of 0.07 eV (from 1.43 eV of the 1^-82^ state of the 
cation to 1.50 eV of the 1^-Bi„ state of the anion, data in organic solid [18a], 
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sec Tabic 4) and 0.11 eV (from 1.31 eV of the 1^A„ state of the cation to 
1.42 cV 1^53^ state of the anion, data in Nc-matrix [29c], see Table 5). when 
going from the cation to the anion of tetracene and pentacene, respectively, 
compare favorably with our computed blue-shifts of 0.07, from 1.70 to 1.77 eV, 
and 0.10 eV, from 1.50 to 1.60 eV (0.05 and 0.08 in Ref. [29^). However, this 
is not the case for the lowest-lying x bands, for which our theoretical pre- 
dictions seem to present a mismatch. In the case of tetracene and pentacene, 
e.g., the experimental blue-shifts of 0.04 eV (from 1.66 eV of the l^B^g state 
to 1.7 eV of the I'^A^ state. Table 4), and 0.10 eV (from 1.27 eV of the I'^B^^ 
state to 1.37 eV of the 1^-623 state. Table 5), are predicted in this study to be 
red-shifted by -0.10, from 1.70 to 1.60 eV, and -0.09 eV, from 1.25 to 1.16 eV, 
respectively (-0.06 and -0.05 in Ref. [29f]). 

We confirm the system-size-dependent errors found for the short-polarised 
transitions {p bands) in the neutral systems [26g,26k]. More specifically, from 

a comparison with the available experimental data reported in Tables 2-6 we 
find that while the relative error in the position of the a band is always of 
the order of 10%, in the case of the p band this error increases from 2% for 
naphthalene to 20% for hexacene (see the triangles on the top panel of Fig. 11). 
On the other hand, we find the position of the /3 band to be reproduced with 
a good precision with relative errors of, at most, 4%. In the case of the radical 
anions and cations, we could not find clear trends as to the performances of 
our TD-DFT/B3LYP results. On the average, the relative errors are larger 
for the long-polarised bands compared to the short-polarised ones for both 
anions (12% vs. 9%) and cations (16% vs. 7%). 

The oscillator strengths / we obtain for the oligoacenes arc also found to 
display systematic changes as a function of molecular size (Fig. 10). Unlike 
the neutral molecules, in all charged species the lowest parallel (y) transi- 
tions have larger /-values compared to the corresponding lowest perpendic- 
ular (x) ones. The oscillator strengths of the x-polarised transitions display 
small changes with the number of benzene units being always of the order of 
0.01 for cations and anions, and decreasing from 0.06 to 0.03 and from 0.09 to 
0.03 for neutrals and dications, respectively. The oscillator strengths obtained 
for the ^/-polarised transitions are found to increase with the molecular size 
for all the charge states considered. In particular, while the increase for the 
neutral molecules is from 1 ■ 10~^ to 1 ■ 10~^ when going from naphthalene to 
hexacene, for the charge states -1, +1, and +2 the corresponding values go 
from 0.03 to 0.37, 0.05 to 0.32, and from 0.09 to 0.65. Therefore, analogously 
to singly-charged species, doubly-ionised PAHs have strong absorption fea- 
tures in the near-IR, visible, and near-UV spectral ranges, a result that might 
be relevant in the astrophysical context. 

Orbital energy differences are well-defined zeroth-order approximations to 
electronic excitation energies [33j]. As shown in Table 7, the HOMO-LUMO 
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energy gap obtained directly as difference of Kolin-Sliam eigenvalues gives a 
better description of the optical gap, i.e., the position of the p-band, with rel- 
ative errors in the range 1-6%, compared to TD-DFT. As already discussed, 
in this latter case a system-size-dependent error is known to exist [26g], with 
an increase from 2 to 20% in the relative error as found in this study when go- 
ing from naphthalene to hexacene. The ASCF QP-corrected HOMO-LUMO 
gaps of neutral oligoacenes as computed at the B3LYP/6-31+G* level, com- 
pare very well with the DFT-based tight-binding GW calculated ones inde- 
pendently from the size of the system, with discrepancies no larger than 3%. 
These values, however, when compared to experiments appear to be affected 
by the same systematic error found for the TD-DFT results (see the diamonds 
on the top panel of Fig. 11). On the other hand, since these errors cancel each 
other in the evaluation of the exciton binding energy, we obtain an accurate 
estimate for Ebind- Appreciable excitonic effects due to both quantum confine- 
ment and reduction of screening are found, with values ranging from 3.9 eV 
for naphthalene to 2.8 eV for hexacene. 



5 Concluding remcirks 



We presented a systematic theoretical study of the five smallest oligoacenes, 
i.e., naphthalene, anthracene, tetracene, pentacene, and hexacene in the charge 
states most relevant for astrophysical applications, namely -1, 0, and -\-2. 
From the ground-state structural relaxations performed at the B3LYP / 6-31-l-G* 
level we computed the electron affinities, the first and second ionisation en- 
ergies, the quasiparticle correction to the HOMO-LUMO gap of the neutral 
systems, and an estimate of the excitonic effects in this class of compounds. 
Good agreement is found with the available experimental data as well as with 
previous theoretical results. To study the electronic absorption spectra we 
used a compendium of the TD-DFT theoretical framework in both real-time, 
to obtain the whole photo-absorption cross-sections in a single step, and fre- 
quency space, to study general trends as a function of charge-state and molec- 
ular size for the lowest-lying valence tt — > tt* in-plane long-polarised and 
short-polarised electronic transitions. The main step forward achieved in this 
work with respect to previous theoretical analyses lies (i) in the spectral range 
covered, that extends up to the far-UV for both neutral and charged PAHs, 
and (ii) in the first detailed study of doubly-ionised species, largely unexplored 
so far. The interest on PAH dications by the astrophysical community has been 
recently renewed by the proposal that they could be plausible candidates to 
explain a red fiuorescence observed in many interstellar sources [32]. We find 
that dications, like their singly-charged counterparts, display strong electronic 
transitions of tt — > tt* character in the near-IR, visible, and near-UV spec- 
tral ranges. As expected, the broad plasmon-like structure peaking at about 
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17.5 cV is found to be relatively insensitive to the charge-state of the molecule, 
but we interestingly find a systematic decrease with increasing positive charge 
of the absorption cross-section between about 6 eV and about 12 eV. Since the 
latter spectral signature is a general property of all PAHs [23c], a comparison 
with astronomical extinction curves could provide an additional observational 
handle for estimating the average charge state of interstellar PAHs. 
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